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Motivated by the recently discovered microwave-induced "zero-resistance" states in two- 
dimensional electron systems, we study the microwave photoconductivity of a two-dimensional 
electron gas (2DEG) subject to a unidirectional static periodic potential. The combination of 
this potential, the classically strong magnetic field, and the microwave radiation may result in an 
anisotropic negative conductivity of the 2DEG. Similar to the case of a smooth random potential, two 
mechanisms contribute to the negative photoconductivity. The displacement mechanism arises from 
electron transitions due to disorder-assisted microwave absorption and emission. The distribution- 
function mechanism arises from microwave-induced changes in the electron distribution. However, 
the replacement of a smooth random potential by the unidirectional one, leads to different relative 
strengths of the two contributions to the photoconductivity. The distribution function mechanism 
dominates the photoconductivity in the direction of the static potential modulation, while both 
mechanisms contribute equally strongly to the photoconductivity in the perpendicular direction. 
Moreover, the functional dependence of the negative photoconductivity on the microwave frequency 
is different for the two mechanisms, which may help to distinguish between them. In another marked 
difference from the case of smooth disorder, the unidirectionality of the static potential simplifies 
greatly the evaluation of the photoconductivities, which follow directly from Fermi's golden rule. 

PACS numbers: 73.40.-c, 73.50.Pz, 73.43.Qt, 73.50.Fq 



I. INTRODUCTION 



Recent experimentsiSi^*^ on a two-dimensional electron gas (2DEG) in weak magnetic fields under microwave 
irradiation have led to the unexpected discovery of regions in magnetic field where the longitudinal resistance is very 
close to zero. Unlike for quantized Hall states, the Hall resistance remains essentially classical and non-quantized for 
these novel "zero-resistance states." These states occur near magnetic fields where, up to an additive constant, the 
microwave frequency w is an integer multiple of the cyclotron frequency lOc- 

This discovery initiated a flurry of theoretical activity from which the following basic picture emerges. It has been 
argued^'^ that under microwave irradiation, the microscopic diagonal conductivity can become negative. This would 
lead to a (macroscopic) instability towards a current carrying state. Macroscopic resistance measurements on this 
state show zero resistance because current can be made to flow through the sample by a rearrangement of large current 
domains. Different microscopic mechanisms for a negative contribution to the microwave-induced photoconductivity 
have been proposed. One mechanism relies on disorder-assisted absorption and emission of microwaves^iSiiSiii (see 
also Ref . Il^ . Depending on the detuning Atj — uJc~uj, the displacement in real space associated with these processes 
is preferentially in or against the direction of the applied dc electric field. In an alternative mechanism, microwave ab- 
sorption leads to a change in the electron distribution function, which can result in a negative photoconductivity>^*i^ii^ 
Detailed calculations within the self-consistent Born approximation suggest that the latter mechanism is larger by a 
factor Tin/r* where Tin is the inelastic relaxation time and r* denotes the single-particle elastic scattering time in a 



* Permanent address 



2 



magnetic field. 

In the present paper, we study the microwave-induced photoconductivity within a model in which the 2DEG 
is subjected to a unidirectional and static periodic potential. Our motivation for doing so is twofold. First, the 
study of periodically modulated 2DEGs in a perpendicular magnetic field has led to the discovery of a number of 
interesting effects such as transport anisotropics^^ and commensurability effects such as the Weiss oscillations of the 
conductivity»iS*i2'iS'iS In addition, the periodic potential lifts the Landau level (LL) degeneracy. This allows one to 
exploit the familiar relation between momentum transfer and distance in real space in high magnetic fields to compute 
the current by applying Fermi's golden rule. In this way, one finds in the absence of microwaves that scattering from 
the disorder potential U leads to a current 

= £^ E E(^^' - k)el\{n'k'\U\nk)m, ~ Jlk'W^k - (1) 

^ ^ nn' kk' 

for an applied dc electric field in the z-direction. Here, is the Fermi-Dirac distribution function of the Landau 
level states \nk) which remain approximate eigenstates even in the presence of the periodic potential (n is the LL 
index) . The (5- function involves the energies e„fc of these states including the effect of both periodic potential and dc 
electric field. 

It is evident from Eq. ^ that the microwaves will affect the current in two ways, (i) The joint effect of disorder and 
microwaves can give additional contributions to the transition matrix elements. This is the origin of the displacement 
photocurrent which relies on the displacements in real space associated with disorder-assisted absorption and emission 
of microwaves. In more conventional terms, this contribution can be associated with the effect of the microwaves on 
the collision integral in a kinetic equation, (ii) The microwaves will also result in a redistribution of electrons, changing 
the electron distribution function /„fe away from its equilibrium form /^j,. This distribution-function contribution to 
the photocurrent will be important if inelastic relaxation is sufficiently slow. Our model allows us to compute the 
various contributions to the photocurrent straight-forwardly within Fermi's golden rule. 

For the parallel photocurrent (i.e., parallel to the wavevector of the static periodic modulation) we find that the 
distribution- function mechanism gives a larger contribution than the displacement mechanism, by a factor Tjn/r*. In 
addition, we find in this case that our results, with suitable identifications, are parametrically consistent with earlier 
results for disorder broadened Landau levels in the self-consistent Born approximation^ii^ By contrast, we find a 
strong enhancement of the displacement mechanism for the perpendicular photocurrent so that in this case, both 
contributions are of the same order. 

This paper is organized as follows. In Sec.|nl we introduce the model and collect the relevant background material. 
In Sec, mil we compute the dark conductivity. The displacement mechanism for the photocurrent is discussed in Sec. 
IIVI while the distribution-function mechanism is worked out in Sec.0 Sec.^also contains a discussion of the Weiss 
oscillations of the photocurrent. The polarization dependence is considered in Sec. IVII We summarize in Sec. IVIII 
Some technical details are given in a number of appendices. In the remainder of this paper, we set h = 1. 

II. THE MODEL 
A. Basics 

In this section, we specify our model and review some relevant background material. We consider a two-dimensional 
electron gas (2DEG) subject to a perpendicular magnetic field B and a unidirectional static periodic potential 

V{r) = Vcos{Qx) (2) 

with period a — 2'n/Q. The periodic potential which hfts the Landau level degeneracy, is assumed to be stronger than 
the residual disorder potential U{v). The disorder potential is characterized by zero average and variance 

{U{y)U{v'))=W{v-v'). (3) 

For white-noise disorder, W{r — r') = ^J"^^ (5(r — r') with Fourier transform M^(q) = l/27ri^T. Here, v denotes the 
density of states at the Fermi energy in zero magnetic field and r is the zero-field elastic scattering time. For smooth 
disorder potentials, the correlator W{y) falls off isotropically on the scale of the correlation length ^ of the disorder 
potential ^ Af for smooth disorder; Xp denotes the zero- field Fermi wavelength). We also note that the impurity 
average of the disorder matrix element between oscillator states \nk) for electrons in a magnetic field in the Landau 
gauge is 

\{nk'\U\nk)\^ = j ^5,^,,,_,e^[L„(^)]^iy(q). (4) 
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Here, n denotes the LL quantum number and k the momentum in the y direction. Ln{x) denotes the Laguerre 
polynomial and is ~ {h/eBY^^ the magnetic length. 

In this paper, we focus on the regime of high Landau levels so that Xp ^ ^ Rc- (Here, Rc denotes the cyclotron 
radius.) We assume that the period a of the modulation satisfies the condition 

AF<a<i?c- (5) 

This is essentially a technical condition, which simplifies some of the calculations. For smooth disorder, we assume, 
in addition, that the correlation length ^ of the disorder potential satisfies the inequality 

AF<C<4/a- (6) 

Here, the first inequality reflects the fact that the disorder is smooth, while the second inequality ensures that the 
typical jump in real space of length £%/^ associated with a disorder scattering event is large compared to the period 
a of the periodic modulation. 

The 2DEG is irradiated by microwaves described by the electric potential 

^(r, t) = -|r(E*e*'^* + Ee"*'^*) = + cj^^e"^', (7) 

where (/>+ = [4>-]* = — eEr/2. The complex vector E contains both strength and polarization of the microwaves. 

In the absence of disorder and microwaves, and for sufficiently weak periodic potential, the single-particle spectrum 
of the electrons can be obtained by treating the periodic potential perturbatively. Starting with the oscillator states 
|nfc), one obtains 

e"^,^u;,{n+^) + VnCOs{Qkil). (8) 

The amplitude Vn is given by 

K = Ve-'^''-/'L^iQ^£l/2). (9) 

In the limit of high Landau levels, Vn can be approximated as Vn ~ VJo{QRc) and thus exhibits slow oscillations 
with period kpa 3> f as a function of LL index n. (The LL index n enters via the cyclotron radius.) This also implies 
oscillations of Vi as function of the magnetic field. It is these oscillations of Vi which arc responsible for the Weiss 
oscillations'''' of the conductivity. 

If in addition, a dc electric field i?dc is applied in the a:;-direction, the eigenenergies take the form 

e„fc ~ fOfc - eE^Ml- (10) 

It is useful to define the density of states (DOS) of a periodic potential broadenend LL by 

^.*(6)=Z.*^*(6) (11) 

with the density of states at the band center 

1 1 



27r£| nVn 



(12) 



and the normalized density of states 



r(e) - . (13) 

v/1 - [{e - En)/VnY 

Here, n and e satisiy |e — i?„| < Vn (with the LL energy En = uJc{n-\- 1/2) ). Note that the DOS can also be expressed 
as V* u{ijJclVn), reflecting the increased density of states due to the Landau quantization. 



B. Kinetic equation 



We now turn to setting up the kinetic equation for the non-equilibrium electronic distribution function fnk which 
describes the occupation of the LL oscillator eigenstates \nk). These occupations change due to disorder scattering. 
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disorder-assisted microwave absorption and emission, as well as inelastic relaxation which we include within the 
relaxation-time approximation. Note that in principle this distribution function depends also on the spatial coordinate 
y. However, it will be sufficient throughout this work to consider distribution functions which are uniform in the y 
direction. (The dependence on x, on the other hand, is included, as the momentum k also plays the role of a position 
in the x-direction.) 

If the dc electric field points in the x-dircction, the kinetic equation takes the form 



dfnk _ ( dfnk \ , f dfnk \ /nk ~ /nk 



dt V dt 1,..^ \ dt 



(14) 



In principle, there should also be a term which describes the drift in the y-direction induced by the dc electric field. 
However, this term has no consequences when considering distribution functions which are independent of y. In the last 
term on the right-hand side, f^j. denotes the equilibrium Fermi-Dirac distribution and Tin denotes a phenomenological 
inelastic relaxation rate. The collision integral for disorder scattering is explicitly given by 

dfr 



f)f I - A^^^\i^'k'\U\nk)\'[U,k' - fnk]5ienk-en'k'). (15) 
The collision integral for disorder-assisted microwave absorption and emission is 

The precise nature of the operator T^. will be given in Eq. 142|) below. Note that these collision integrals involve the 
electron energies including the effects of the dc electric field. 

If the dc electric field points along the y-direction, we can no longer include it in the eigenenergies. Instead, it 
enters the kinetic equation through an additional term describing the associated drift in the a;-direction, 

dfnk _ dfnk , f dfnk \ , f dfnk \ /nk ^ /nk 



The collision integrals are given by the expressions in Eqs. H15|l and (|16() with the energies in the (5-functions taken in 
the absence of the dc electric field. 

We close this section with a calculation of the elastic scattering rate l/r* in high magnetic fields. The motivation for 
doing this is twofold. First, r* is a natural parameter in terms of which to write our final results for the conductivity. 
On a more technical note, computing r* gives us the opportunity to introduce a convenient way of dealing with 
integrals involving Laguerre polynomials, which will be used repeatedly throughout this paper. From the collision 
integral for elastic disorder scattering, we obtain 



J— ^y^27r\{n'k'\U\nk)\^Sie- Cn'k') (18) 



with e = e„fe. Noting that n ~ n' and inserting the expression Q for the matrix element, we obtain 

^ 2^ / ^e^[LnC^)YW{ci)5{e-enk+,,). (19) 



T*(e) J {2n 

The Laguerre-polynomial factor arising from the matrix elements of the disorder potential decays as a function of q£^ 
on the scale of the cyclotron radius Rc, in addition to fast oscillations on the scale of the zero-field Fermi wavelength 
Xp- On the other hand, the argument of the 5- function changes with on the scale of the period a of the periodic 
potential. Thus, for white-noise disorder and in the limit Xp <^ a ^ Rc, we can average the 5-function separately 
over Qy. Using the identity 

{die- e'nk'))k' ^ 27:11,.* (e) (20) 
and performing the remaining integral over the Laguerre polynomial. 



(27r)2 ^ 2 27r£2j' 
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we find the result 

1 _ 1 v*{e) 
r*(e) T 

In the following, we will also use the notation r* = r*(e = En), i.e. 



(22) 



r*=r^. (23) 

This result reflects the increased density of final states in the limit of well-separated Landau levels. 

For smooth disorder, we need to distinguish between the single-particle scattering time and the transport scattering 
time. Their zero- field values Tg and Ttr are related to the finite field values T*(e) and T^*(e) in analogy to Eq. I|22|l . 
i.e., r*(e) = Tsv/v*{e) and t^* (e) = rtri^/i^*(e). Some details of the calculation are given in Add. 1^ 

III. DARK CONDUCTIVITY 
A. Conductivity (Jj:^ along the modulation direction 

In this section, we compute the dark conductivity, i.e., the conductivity in the absence of microwaves. We start 
with the situation in which the dc electric field is applied in the x-direction, i.e. parallel to the wavevector of the static 
periodic modulation. We assume that the dc electric field is sufficiently weak so that heating effects can be ignored. 
In this case, the distribution function remains in equilibrium, f^k — fnk^ ^^'^ system responds to the dc electric 
field with a current in the x-direction. 

This current can be expressed by counting the number of disorder scattering events that take an electron from 
a state fc, localized in the x-direction at k£g to the left of an imaginary line xq parallel to the y-axis, to a state 
A;', localized at fc'£^ to the right of this imaginary line, and vice versa. Due to current conservation, the current is 
independent of the particular choice of xq and it turns out to be useful to average over all possible xq. This results 
in the expression 

for the current in the a;-direction. Performing the integral over xo gives 

= 7^EE(^' - k)ll\{n'k'\U\nk)m^ - f°,,,]S{enk - (25) 

^ nrv kk' 

Inserting the explicit expression Q for the disorder-averaged matrix element and performing the sum over k\ one 
obtains 



^ nk 



E / -0j2ly^Be^[Lni^)rWiq)[f^, - fnk+JS{el, - 6°,+,^ + ei?dcg,4)- (26) 



Expanding to linear order in the dc electric field, one obtains for the conductivity 



^ nk 

For white- noise disorder and Xp ^ a <ti Rc, the Laguerre-polynomial integral can be computed in analogy with the 
evaluation of Eq. 1)19(1 above. Using 



this yields 



T2(%4)^e^[i«(^r = ^, (28) 



^^"^ ' t.I-'tS')'-'^*^^'^- (29) 



L^Ly 2-KVT ^ \ de^^ 



ik 
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Expressing the sum by an integral involving the density of states, we can cast this result in the final form 

df{e) 



cTxx = / rfe ( ^— ) axxie) (30) 



(31) 



in terms of 

V2rt*(e) 

This equation is written such that it includes both types of disorder. For white- noise disorder r^* ~ r* ~ r*, while 
for smooth disorder r^* ^ r*. The derivation of the result for smooth disorder is sketched in Add. IXI 

This result for the dark conductivity can be interpreted as follows. The bare rate for disorder scattering is l/r*, 
where each scattering event is associated with a momentum transfer 1/^. This momentum transfer translates into 
a jump of magnitude in real space so that the electron diffuses in the z-direction with a diffusion constant 

Dxx ^ {^'b / 0'^ / 'Ts ■ Alternatively, this diffusion constant can be written in terms of the transport time as Dxx = 
RI/2t^j. (using that t^^/t* ^ {kp^,)'^)- By the Einstein relation, this diffusion constant translates into the conductivity 
given in Eq. The conductivity can also be expressed in terms of the zero-B conductivity axx{B = 0) as 

CTxx — <^xx{B = Q) / {lUcT^^Y . We also note that axx ~ 1/^™ so that the oscillations of V with magnetic field B (cf. 
Eq. lO above) lead to Weiss oscillations of the conductivity, in agreement with previous resultsJ^ 

The energy integral in Eq. H3U|I is formally logarithmically divergent due to the square-root singularity of the density 
of states v*{() at the band edge. This singularity is cut off by smearing of the band edge by disorder or by the applied 
dc electric field, when the latter is kept beyond linear order. 



B. Conductivity (jyy perpendicular to the modulation direction 

An applied dc electric field in the y-direction leads to a non-equilibrium distribution function /„fc due to the drift 
term in the kinetic equation H17|) . In the absence of microwaves, linearizing the stationary kinetic equation in the 
applied dc electric field yields 

ei?dc^-2vr j ^e^[LX^)fW{ci)[5U+,,-5U]5{el,^el^^J. (32) 

Here, 5fnk = fnk — Ink denotes the deviation from the equilibrium distribution function, and we have neglected 
inelastic processes relative to elastic disorder scattering. Due to the periodicity in the x-direction, 5 fnk = Sf^k+a/e'^g ■ 
Moreover, if k and k + qy are two momenta with the same energy 6,^^^,, but opposite signs of the derivative de^i^/dk, 
then Sfnk = —Sfnk+qy- Using that for white-noise disorder and Xp ^ a ^ Rc, we can split the q-integration as in 
the evaluation of ifT^ and obtain 

6U = -eE^cr*{el,)^. (33) 
In terms of the distribution function, the current in the y-direction is given by 

Cfc„i, 



3y 



e 



^x-^y 



ik 



Inserting the expression for the distribution function gives for the conductivity 



'yy 



de[-^4r^]cTyy{e) (36) 



Expressing the sum over nk as an energy integral involving the DOS, we obtain 

n 

with 

ayy{e) = emvy{e)rr:{e))i^*{e). (37) 
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As above for a^xi this result is written such that it includes both the case of white-noise and of smooth disorder. The 
derivation for the case of smooth disorder is sketched in App. ^ We have defined the drift velocity 



dk 



^ (38) 



in the y-direction, induced by the periodic modulation. 

The result H37|l can be interpreted as follows. With respect to the motion in the y-direction, a partially filled LL 
consists effectively of a set of two "internal edge channels" parallel to the y-axis per period a. Neighboring channels 
flow in opposite directions so that disorder scattering randomizes the direction of the motion in the y-dircction after 
time T*. The factor Dyy — VyT* can thus be interpreted as the diffusion constant of the resulting diffusion process. 
We note that unlike (7xx{^)i the conductivity a-yy{e) remains finite at the band edge. The anisotropy (Tyy/cxx of the 
dark conductivity is thus of order {vyT* / RcY {kp^Y , where both factors are larger than unity. The dark conductivity 
Gyy depends on the modulation- induced LL broadening as Uxx K?' that the Weiss oscillations in Gyy are phase 
shifted by tt relative to the oscillations in axx^ in agreement with standard resultsi^ 

The dc electric field also leads to heating of the electron system. The characteristic field E*^^ where this becomes 
relevant, can be estimated as follows. The dc electric field causes a drift in the x-direction with drift velocity (Edc/B), 
changing the potential energy of the electron by {V/a){Edc/B)T* . This gives rise to a diffusion constant in energy of 
£>£ ~ {y / af'{Edc/ B^T* . Heating can be neglected as long as the typical energy change (ZJcTin)^/^ is small compared 
to V . This gives the condition 

£;dc « El, = (39) 

for the dc electric field, ft is only for these electric fields that the result H37() is valid. A more formal derivation of this 
result is given in App. El 

For larger dc electric fields Edc ^ E^,, the effect of heating needs to be taken into account. Following the arguments 
given in App. ^ one expects that the conductivity is suppressed by heating effects and behaves in magnitude as 

^^'e^DyyU*. (40) 
The reason for this suppression is that heating reduces the fc-dependence of the distribution function. 



IV. DISPLACEMENT PHOTOCURRENT 



A. t-matrix elements 



The microwaves lead to additional contributions to the transition matrix element between LL oscillator states 
which enters into the current expression Direct microwave absorption or emission does not contribute to the 
current, because the microwaves do not transfer momentum to the electrons so that such processes are not associated 
with displacements in real space. In addition, such processes occur only for lj = lJc- On the other hand, disorder- 
assisted microwave absorption and emission is associated with displacements in real space of the order of Rc {£% /£, 
for smooth disorder). This process is allowed for microwave frequencies away from lUc- In this section, we compute 
the contribution of this displacement mechanism to the photoconductivity within our model. 

The transition rate between LL oscillator states involves the i-matrix 

T= ([/-K ,/)) + ([/ + (/))Go ([/ + (/)) + ..., (41) 

where Go denotes the retarded Green function of the unperturbed system. The dark conductivity, computed in the 
previous section, follows in the approximation T ~ [/. Disorder-assisted microwave absorption T+ and emission T_ 
is given by 

T± ^ [UGo(^± + (b±GoU]. (42) 

Note that T+ and r_ contribute incoherently. Assuming that lo couples only neighboring LLs and that the microwaves 
are linearly polarized in the x-direction, the corresponding matrix elements between LL oscillator states are 

{n±lk'\T±\nk) ^ ±(^^)[{n±lk'\U\n±lk)- {nk'\U\nk)], (43) 
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where we used 

Go,n±ife(e„fc ±Lu) — j = T-7— 

ink — Cnztlfc 



Go.nk'ienk) = =±x-- (44) 

Using the disorder matrix elements and neglecting corrections of order 1/n, one finds 

|(n±lfc'|r±|nfc)p^ (^yi _05^^,„_,e-^[W(^)-L„(^)mq). (45) 

Thus in this approximation, the matrix elements are identical for absorption and emission and depend only on the 
absolute value oi k — k' . It is worthwhile to point out that the difference between the two Laguerre polynomials 
reflects the fact that disorder-assisted microwave absorption involves a coherent sum of two processes: In one process, 
a microwave photon is first absorbed, resulting in a transition from the nth to the {n + l)th LL, with disorder 
subsequently inducing a transition between states in the {n + l)th LL. In the second process, disorder first leads to a 
transition between states in the nth LL with a subsequent absorption of a microwave photon. The divergence of the 
matrix element for Aw ^ is an artefact of low-order perturbation theory in the disorder potential U . In a more 
accurate treatment, this divergence would be removed by disorder broadening. 

B. Displacement photocurrent jx along the modulation direction 

In this section, we compute the displacement contribution to the photocurrent in the modulation direction. The 
current in the .T-direction can now be computed in terms of Fermi's golden rule in the same manner as for the dark 
current in Sec. Illll In this way, one obtains 

^TJ-T. E(fc' - + lA:'|r+|nfc) - - 6„+i,, + w). (46) 

X y 



n kk' 



Here, we assume that the microwaves frequency w > is such that it couples neighboring LLs. 

In order not to complicate the calculations unnecessarily, we will consider temperatures T ^ V. In this regime, the 
temperature smearing is over an energy range large compared to the LL width and the distribution function depends 
only on the LL index n, but not on the momentum k. Inserting the explicit expression (|45|1 for the i-matrix element, 
we obtain 

■photoi 27r£|e f eER^ Y^ 



X 

k 



E / (^9.e-^[L„+i(^)-L„(i|^rH^(q)<5(e„+i, -a;). (47) 



The fc-summation gives 



E - - -) - 2<2\[3in^^-(A|)2]i/2' (48) 



where we introduced AtD = Aw — eE^>^c%^%- Thus, we obtain for the current 



photoi _ ^ I t^^J^c \ \^rr _ r 1 



2Tr£l V 4Aw 



d^g 2-i!|ir, (I'el. q'il , Wjq) 



[2^^'^^^ 2 ' 2 v:„[sin2%-(^)2]i/2 

where the integral is only over the region where the square root in the denominator is real. It is useful to interpret the 
various factors in this expression. It consists of a charge density per LL e/27r£^, and a rate (per LL) for jumps in the 
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FIG. 1: The functions Ai{Aio/2VN) (full line) and Bi{Alo/2Vn) / yn{VN / A) (dashed line), describing the dependence of the 
parallel photoconductivity on the microwave frequency, cf. Eqs. I|52^ and 1)68^ . 



x-direction with lengths between qyl'^ and {qy + dqy)Pg, multiplied by the jump lengths qyi% Finally, the expression 
is integrated over all jump lengths qy and summed over all LLs. 

The sum over LLs n is trivial and for white-noise disorder, the integral over q can again be decoupled for Af <C 
a <ti Re- Expanding to linear order in the dc electric field and using the integral 



(2^)2 ^""■"^^ - L—T.V 2 ' ""^ 2 " TT 

we obtain the linear-response conductivity 



(7: 



photo I 



XX 



eERc 



2 



A,iAu;/2V^) (51) 



with the function (cf. Fig. ^ 



3 d 



{x) = - — —K ^Jl-xA , (52) 



expressed in terms of the complete elliptic function K. Note that the first factor in Eq. (|51|l is just the dark conductivity 
axx [cf. Eq. H31|l ] and that the result includes the case of smooth disorder. (Strictly speaking, the result for smooth 
disorder is valid only up to a numerical prefactor that depends on the precise nature of the smooth disorder potential, 
cf. Add. 1X1 for details.) The behavior of this displacement photocurrent for Aw ^ 2Vn and Acj ~ 2Vn follows from 
the asymptotic expressions 

i^(v/l-(Aw/21V)2)c.| ^{1 + ^) «^i_(A^)2«i • (53) 

This implies that Ai{Au!/2Vn) remains finite for |Aa)/2V;v| ~^ 1 (cf. Fig. [21 and is proportional to 1/Alu for small 
|Ati;/2V;v|- The sign of the displacement photocurrent is given by sgn(a;c — to), similar to previous work on disorder- 
broadenend LLs&ii 

The magnitude of the displacement contribution H51|l to the photoconductivity can be understood as follows. The 
bare rate of disorder-induced microwave absorption is {l/T*){eERc/ Auj)'^ , where the second factor is the dipole 
coupling of the microwave field for Landau states divided by the relevant energy denominator of the intermediate 
state. Each of these scattering events is associated with a jump in real space of the order of resulting in an 

effective diffusion constant Dxxi^ERc/ AloY . An additional factor (r*/r*j,) arises because of the partially destructive 
interference of the two contributions which were discussed below Eq. (|45() . This interference leads to the difference 
of Laguerre polynomials in Eq. 1)49(1 which introduces an additional factor [q/kpY into the integral. This factor is of 
order l/ikpO^ ~ (r,Vr*). 
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FIG. 2: Illustration of disorder-assisted microwave absorption for Atj ~ 2Vn- In the a;-direction, the Landau levels are 
modulated by the periodic potential and tilted by the dc field. The square-root singularity of the Landau level DOS 11311 at 
the band edge leads to the singular behavior of the photocurrent for these detunings Alj, cf. Fig. Q 



C. Displacement photocurrent jy perpendicular to the modulation direction 

The current in the y-direction can also be computed in a semiclassical approach. If a dc electric field is applied in 
the y-direction, the equipotential lines of energy £ are "meander" lines defined by £ — Vn cos{Qx) — eE^cV- This gives 

y= ^[VnCosiQx)-£] (54) 

with an average y- value y = —£/eEdc- Quantum-mechanically, we can think of these equipotential lines as states. 
This relation as well as the calculation sketched in this section is worked out more formally in App. |0 Scattering 
between such meander states that differ in energy by Aw = ujc — uj (ignoring the LL energy) involves jumping a 
distance Ato/eEdc in the y-direction. Thus, we can again compute the current by Fermi's golden rule. It is important 
to observe that the direction of the jumps is fixed by the sign of the energy difference. Below, we will comment on 
the limits of validity of this approach. 

The current expression involves the rate of jumps. If Edc is sufficiently weak, the amplitude of the meander line is 
very large compared to the scale Rc over which jumps occur. On the scale of the jump, the meander lines are therefore 
essentially indistinguishable from the equipotential lines in the absence of the dc-field. This allows us to employ the 
rate of jumps which we obtained from the calculation of the current in a;-direction. (We have to set E^c = in the 
formulas obtained there, cf. Eq. (|49|l l. In this way, we obtain for the displacement photocurrent in the y direction 

■photoi „ e f eERc VM-^ 



As mentioned above, we derive this expression more formally in App. |21 Computing the integral for white-noise 
disorder, we obtain the result 

= '^1 ( ^V^^(v/l-(A-/2W), (56) 



•'y {27T)WNi%Edc V 4Aa; J 2'kvt 

valid for Xp <ti a <ti and T 3> Vpf. Note that the current is not linear in the applied dc electric field but rather 
diverges as as l/E^c- This anomalous behavior is associated with the fact that the length of the jumps diverges with 
decreasing dc electric field and that the direction of all jumps is the same, fixed by the sign of Ato. Rewriting the 
result in terms of a conductivity, we obtain 

= l^'D,.y] [^^^^J (^)^,(A„/2V„). (57) 

Note that the first factor is just the dark conductivity ayy. The derivation of the result for smooth disorder is sketched 
in Add. IXI The function ^2(0;) is defined by 



A2{x) = 2xK{^/l-x^) 



(58) 
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FIG. 3: The functions y42(Att;/2VAr) (full line) and i32(Att;/2VAr) (dashed line) describing the dependence of the perpendicular 
photoconductivity on the microwave frequency, cf. Eqs. I58|l and Ij71|l . 



and plotted in Fig. The sign of the photocurrent is given by sgn(tt'c — ^\ as in the case of a?^^°^ . 

The magnitude of 0-^^°*° can be understood as follows. Since all jumps are in the same direction, we es- 
timate the current density directly. Effectively one LL contributes so that the relevant density of electrons is 
l/27r£g. The step length is Aw/ei?dc and the rate of jumps is given by (l/r*)(ei?i?c/Ai:t;)^(T*/T^*) where the 
first two factors are the bare rate and the last factor again reflects the partially destructive interference between 
the two contributions to disorder-assisted microwave absorption. Thus, we find a current of order ]^°'^°^ ^ 
{e/2'Kl%){Auj/eEAc){l/T*){eERc/ Aujf{T* /t^^), in agreement with Eq. (|^. 

The limits of validity of this result are most naturally discussed in terms of a semiclassical picture. Semiclassically, 
the individual scattering events such as disorder-assisted microwave absorption leave the j/-coordinate of the electron 
essentially unchanged (to an accuracy of Rc)- The full jump by Auj/eE^c is realized only if the electron remains in the 
meander state it scattered into for sufficiently long times to explore its entire y-range. Under the condition t* <^ Tin, 
the electron will diffuse on the meander line before equilibrating by inelastic processes. The typical diffusion distance 
y/DyyTin in the y-direction should be larger than the amplitude V/v/e-Bdc of the meander line. Thus, we find that the 
condition for the validity of the expression (|57|l is 

i?dc » = (59) 
27rVr*ri„ 

Note that this is just the opposite of the range of validity of the dark conductivity ayy computed in Eq. H37|l . 

For smaller dc electric fields, the jumps are no longer all in the same direction and we can estimate the displacement 
photoconductivity as follows. The disorder-induced microwave absorption excites the electrons to a meander (equipo- 
tential) fine which is shifted in the y-direction by Auj/eEdc relative to the initial state. For definiteness, assume that 
this shift is in the positive y-direction. Since quasiclassically, the jump itself leaves the y-coordinate of the electrons 
unchanged (to an accuracy of Rc), the electrons will initially populate only those parts of the excited meander line, 
which are at least a distance Aui/eE^c from its top (in the y-direction). After the excitation, the electrons begin to 
diffuse on the equipotential line due to disorder scattering, typically a distance \jDyyT\i^ before they relax back. Thus, 
after time Tin, the population of the excited meander line will extend further in the positive y-direction by a distance 

DyyT\^^ and the average positive drift per electron is [\/ DyyTi^/ {Vn / eE(ic)]\/ DyyTi^. These arguments lead to a 
displacement photoconductivity of 

(60) 

valid for E^^ *C ^'dc- ^^ote that this result matches with Eq. H57|) for E^^ — E^^. It is interesting to note that even 
the linear-response displacement photoconductivity involves the inelastic relaxation time Tin- 
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V. THE EFFECT OF A NON-EQUILIBRIUM ELECTRON DISTRIBUTION ON THE 

PHOTOCURRENT 

A. Distribution function 

For non-zero inelastic relaxation time Tin, the microwave irradiation changes the electron distribution function, 
away from the equilibrium Fermi-Dirac distribution. In this section, we consider the contribution to the photocurrcnt 
arising from this change in fnk- 

The microwave-induced change in the distribution function can be computed from the stationary kinetic equation 
in the absence of the dc electric field. Elastic disorder scattering contributes only when states of the same energy have 
different occupations. Since this is not the case for iJdc — 0, we can ignore elastic disorder scattering when computing 
the microwave- induced change in fnk- Thus, the kinetic equation reduces to a balance between the microwave-induced 
collision integral and inelastic relaxation which yields to linear order in the microwave intensity 

6 fnk = fnk - fL - n„ E E 2vr|(n'fc'|T,|nfc)n/^,,, - /0,]5(e^, - e^,,, -I- au). (61) 

n'k' cr— ± 

As before, we restrict attention to temperatures T ^ ^ so that f^^ — fn^ independent of k. Inserting the explicit 
expression H45f) for the t-matrix element for white-noise disorder and using the decoupling of the q-integration for 
Ai? <C a ^ i?c, we obtain for the change in the distribution function 

Note that k enters this expression only through e° ^ . Strictly speaking, this expression breaks down when e° — aAuj 
approaches the band edge. In this limit, it is no longer sufficient to treat the microwave field to linear order in the 
intensity. Effectively, this divergence is cut off when Sfnk becomes of order unity, i.e., for distances Ae from the band 
edge satisfying 



eER, 



(63) 



Here, we used that the DOS has a square-root divergence at the band edge. We also note that our linear approximation 
in the microwave intensity breaks down completely beyond microwave intensities given by (eERc/ Auj)'^{Tia/T^j.) 1. 
In the estimate 

Sfnk ^ (l/T;)(ei?i?,/Ac.)2(r;/T*)Ti„ (64) 

for the magnitude of Sfnk, the first three factors are the rate for disorder-assisted microwave absorption and emission. 
The last factor Tin represents the time interval during which electrons are excited. The expression (|62|l will form the 
basis of our calculation of the distribution-function mechanism for the photocurrcnt to which we now turn. 

B. Distribution-function contribution to the photoconductivity along the modulation direction 

Going through the same steps as in the derivation of the dark current in the x-direction, one finds that Eq. H29|) 
for axx remains valid even for a non-equilibrium distribution function, as long as it depends on k through e^^^, only. 
Thus, we find for the distribution-function contribution to the photoconductivity 



photon ^ ^2.£|-^ V (-^^ ) .*{el). (65) 



Inserting 5 fnk from Eq. H62(l . performing the sum over the LL index n, and rewriting the sum over k as an energy 
integral involving the density of states, we obtain our result 

^photon ^ 4 f^^V [<^'DxxU*] B,iAu;/2V^) (66) 



4Aw 
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for the distribution-function contribution to the photoconductivity. The function 



Si(Aw/2FAr) = --— / dev*{e+\^uj\)[v*{e)Y. (67) 



Vn 



is of order unity for Alj ~ Vat. This result shows that for the parallel photocurrent, the distribution- function 
contribution is larger than the displacement contribution H51() by a large parameter Tjn/r*. This result is consistent 
with earlier results for disorder-broadened Landau levels»ii*i^ 

The integral H67|l has a logarithmic singularity at the lower limit, which has the same origin as the divergence as the 
dark conductivity Uxx in Eq- (|31|) . Thus, the singularity is cut off in the same way as for the dark conductivity. (For 
small Aw, one may seemingly have a more serious singularity. However, we should remember that the DOS arising 
from (5/ never really becomes singular, cf. the discussion above around Eq. (|63|l .^ Since the logarithmic singularity 
dominates the integral 1)67(1 . we can replace e by the lower limit in the DOS v*{e + |Aa;|). In this way, we obtain the 
result (cf. Figs.QJ 

1 l-2|a;| Vn 

Here, A denotes an effective broadening in energy of the band edge, either due to disorder or a finite dc electric 
field, which cuts off the logarithmic singularity. Interestingly, this shows that the sign of the photocurrent is given 
by sgn([Jc ~ for |Alj| < Vn only. For lAwj — Vn, we find additional sign changes which are associated with the 
singular nature of the DOS at the Landau level edge, cf. Eq. (fT^ . 



C. Distribution function contribution to the photocurrent perpendicular to the modulation direction 



To compute the distribution- function contribution to the photoconductivity Uyy, we note that Eq. ((35|l remains 
valid for non-equilibrium distribution functions which depend on k through e^j,, only. Thus, our starting point is 



KkS" 1 95 Ink 



LxLyj^y dk J ^*(e°fc) denk ' 
Inserting Sfnk from Eq. ((62|l and rewriting the sum over nk as an integral, we obtain the final result 



vv 



Tin \ / eERc 



photon ] ^e^Dyyiy*]B2iALu/2VN). 



The expression (|70f) is given in terms of the function 

g rVN 



B2{Auj/2Vi 



de- 



dAu ./_y„+iA.r^[i>*(eF 



i?*{e-\Au;\). 



This integral is elementary and we obtain (cf. Fig. 



B2 (x) = 



\x\ (^arcsin(l - 2 |x|) + ^) - 8 



sgnx. 



Asymptotically, this function behaves as 



B2ix) 



4v/r 



a;|sgnx 



Ixisgnx 



1- kl < 1. 



(69) 



(70) 



(71) 



(72) 



(73) 



Interestingly, B2{x) has different signs in these two limits, implying that the function -62(2;) must have a node between 
the arguments a; = and x = 1. This node is approximately at \ Aljj\ k, 2Vn /t^"^ ~ Q.2V. As a result, the distribution- 
fimction contribution to the photoconductivity crP^°*°^^ has the same sign as 0"^^°*°^ only in the range \Auo\ > 0.2V;v- 
Again, we associate this behavior with the anomalous behavior of the DOS. 

The magnitude is, apart from a function of Aw/Vat, of order ayySf. Remarkably, in this case the magnitude is of 
the same order as the displacement contribution in Eq. (|60|l . The reason for this is that the displacement mechanism 
exhibits a singular, non-Ohmic dependence on the dc electric field which is cut off at small fields by inelastic processes 
only. A more detailed analysis beyond this order of magnitude comparison would require an accurate calculation of 
the displacement contribution to the photocurrent in the linear-response regime. Such a calculation is beyond the 
scope of this paper and the result may in any case be sensitive to details of the model for inelastic relaxation. 
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D. Weiss oscillations of the photocurrent 



The Weiss oscillations arising from the oscillatory behavior of the Vn as function of LL index n or magnetic 
field have two effects on the photocurrent. First, they lead to a modulation of the amplitude of the photocurrent. 
This amplitude modulation is similar to that of the dark conductivity as the photoconductivity is proportional to 
the dark conductivity. A difference may arise from the additional prefactor Tin/rj* entering the photocurrents, cf. 
Eqs. (|66|l and H7U|) . Specifically, if the ineleastic relaxation rate ri„ depends in a different way on the LL DOS 
compared to the transport time Ttr, there may be a distinct difference between the Weiss oscillations in the dark 
and the photoconductivity. Depending on whether the Weiss oscillations in the photoconductivity are more or less 
pronounced than those in the dark conductivity, this may help or impede reaching negative conductivities and hence 
observing the zero-resistance state. 

A second effect of the Weiss oscillations is associated with the modulation of the LL width Vn- The expressions for 
the photocurrent (|66|l and l|7U|) involve a factor which depends on the ratio of the detuning Auj and the LL width Vn 
at the Fermi energy. This implies that in the limit of well-separated LLs considered here, the range in the detuning 
over which there is a significant photocondictivity also oscillates with magnetic field or Fermi energy. 

In the photoconductivity, the l/B-periodic Weiss oscillations are superimposed on the microwave-induced oscilla- 
tions which are also periodic in 1/B. The periods in 1/B of these oscillations are ea/mvp and e/muj, respectively, 
which can be of comparable magnitude. 



VI. POLARIZATION DEPENDENCE OF THE PHOTOCONDUCTIVITY 



In this section, we discuss the dependence of the photoconductivity on the polarization of the microwave field E. 
We begin by calculating the transition matrix elements for a microwave field polarized in the j/-direction. In this case 

4> = - — y(e*'^* + e-"^*) = + ?!)_e'"* . (74) 



The matrix elements of the operator (l)± are given by 

i dk 



n±mk'\cj)±\nk) = - ( l^Sik' - k)) e'^''-'''^"'"^/* 



,„,LU'-^^)^K...^^Ll(^^^]l ,75) 



2 J ' ""^^ V2n 

for m > 0. Using L\{0) = 1 and L}^{0) = n -I- 1, we obtain for the transition matrix element T± the large-n result 

(n ± lk'\T±\nk) = -i (^^^ j ± lk'\U\n ± Ik) - {nk'\U\nk)] (76) 

valid for Alu <ti u)c- This matrix element differs from the corresponding matrix element for polarization in the x- 
direction by a multiplicative prefactor of unit modulus. As the photoconductivity depends only on the modulus of 
the transition matrix elements, we find that the photoconductivity is the same for microwave fields linearly polarized 
in the x and y-directions. One also readily concludes that more generally, the photoconductivity remains unchanged 
for any linear polarization of the microwaves. 

We now turn to irradiation by circularly polarized microwave fields which are described by 

{{x±iy)e'^-' + {xTiy)e'^') . (77) 



2^/2 



Combining the transition matrix elements for microwaves linearly polarized in the x and y-directions, one finds zero 
photoconductivity for (f)„^ . In this case, the E- vector rotates opposite to the circular cyclotron motion of the electrons 
in the magnetic field. For , both E and the cyclotron motion rotate in the same direction and the photoconductivity 
is double that for linearly polarized microwave fields. 



VII. CONCLUSIONS 



We have computed the microwave-induced photocurrent in the regime of high Landau level filling factors, in a 
model in which the Landau levels are broadened into a band due to a static periodic modulation. Assuming that 
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the static modulation is small compared to the spacing between LLs, the eigenstates are still given by the Landau 
level oscillator states in the Landau gauge. The localization properties of these states allow us to compute the dark 
conductivity as well as the microwave-induced photoconductivity using Fermi's golden rule. The Fermi's golden rule 
expression for the current directly suggests that there are two distinct mechanisms contributing to the photocurrent, 
analogous to previous results^iii^ for disorder-broadenend Landau levels, (i) The displacement mechanism relies on 
the spatial displacements associated with disorder-assisted microwave absorption and emission. This contribution 
can be associated with an additional, microwave-induced contribution to the transition matrix element in the Fermi's 
golden rule expression, (ii) The distribution-function mechanism by contrast, relies on the microwave-induced change 
in the electronic distribution function, again due to disorder-assisted microwave absorption and emission. 

For the photocurrent parallel to the modulation direction, we find that the distribution-function mechanism [cf. Eq. 
(|66|l ] dominates by a large factor Ti^/r* over the displacement contribution [cf. Eq. H51|l ]. in agreement with earlier 
results for disorder broadened Landau levels«ii*i^ The sign of the photocurrent changes with the sign of the detuning 
Auj = — w of the microwaves. For the dominant distribution function mechanism, there are additional sign changes 
associated with the divergence of the density of states at the edge of the Landau level. Remarkably, the situation is 
rather different for the transverse photocurrent perpendicular to the modulation direction. In this case, we find that 
the displacement mechanism is in some sense singular with the result that both contributions to the photocurrent 
[cf. Eqs. (|57|l and (|70() ] are of the same order. We find that our results remain unchanged for any linear microwave 
polarization. For circular polarization, we find a nonzero photoconductivity only when the microwave electric field 
rotates in the same direction as the cyclotron rotation of the electrons in the magnetic field. 
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APPENDIX A: SMOOTH DISORDER POTENTIALS 



As opposed to a white-noise potential, the (zero-magnetic-field) single particle time Ts is different from the transport 
time Ttr for a smooth disorder potential. Specifically, the single-particle time can be expressed in terms of the correlator 
W as 



1 1 

- = 27r^VK(q)(J(efe-efc+,))FS- / dqW{q}, (Al) 

Ts ^ TTVF Jo 

where the average is over the Fermi surface and ek denotes the zero-field dispersion. Likewise, the transport time can 
be expressed as 

1 1 P°° 

- = 27rV(l-cos0q)iy(q)(,5(efc-efe+,))Fs = / dq{q^ /2kl)W{<i), (A2) 

Ttr ^ TTVF Jo 



where 6q denotes the scattering angle. Note that th/ts ^ (fc^?^)^. 

We start by considering the elastic scattering times for smooth disorder. For r*, we need to reconsider the q- 
integration in Eq. ((T^ . 



d^q ,g^^lM2 



Io = J j^e^[L4^)rW{q)S{e'^, - e^.^J. (A3) 

In the limit of large iV, the Laguerre polynomial has oscillations on the g-scale of Xf/(^'b ^'^d falls off on the scale of 
Rc/i%. The correlator falls off on the scale 1/^ and finally, the characteristic scale of the argument of the 5- function 
is a/t\. Thus, unlike for white-noise potential, it is now the correlator M^(q) which cuts off the integral at large q. 
Under the condition Aj? <C a ^ we can still factorize the g-integration as 

h= j ^,e^[L^{^)YW {cornel, - e°,,)>fe'- (A4) 
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Since the integral is now cut off at large q by W, it is sufRcient to employ the semiclassical equation 



e^L^i^) ^ yiJIcosCgi?, - n/A) (A5) 

in the remaining integral. (Strictly speaking, we need a slightly more accurate approximation. However, this changes 
only the argument of the cosine^'' which does not affect the results.) This yields 



dgty(q)((5(e°, - el,,)),, = . (A6) 

ly ZttTs 

The same integral is involved in the computation of the dark conductivity ayy. 
The transport time involves the integral 

l2=J ^^e^[L,.{'^)rW{ci)S{el,-el,^J, (A7) 

where we used that 1 — cos^q ~ q'^/2kp for g ^ kp- An analogous analysis as for Iq above gives the result 

^ (A8) 

1^ 27rrti. 

The same integral appears in the calculation for the dark conductivity axx- 

A similar integral appears in the calculation for the displacement photocurrent f namely 

For smooth disorder, the difference of Laguerre-polynomials is suppressed relative to a single Laguerre-polynomial. 
Using that qRi"''^^^ ~ qRc^'^ + q/kp, we find that the difference effectively introduces a small factor q'^/kp into the 
integrand and thus (cp. the calculation for I2 above) 

Jo^^ — . (AlO) 

For the displacement photocurrent crP^°*°\ we need to consider 

^'^J (0 (9V24)e^[i™+i(^) - LU'^)rW{nm, - (AH) 

Again, the difference of Laguerre polynomials introduces a small factor q'^ /kp into the integrand. The resulting 
integral can no longer be related directly to either Ttr or r^. However, noting that every factor {q/kp)"^ reduces the 
integral by a factor of order l/(fci?^)^, we can estimate 

J, . illM^ 1 . _L4^ (A12) 

The precise numerical prefactor depends on the detailed nature of the smooth potential. 



APPENDIX B: DISTRIBUTION FUNCTION FOR LARGE dc ELECTRIC FIELDS 



For large dc electric fields, E'^jc S> E'^^^ Joule heating effects become important. In this appendix, we study 
the distribution function in this limit. We start by decomposing the distribution function /„fc into symmetric and 
antisymmetric parts (T„j, and a^k under k — > — fc, i.e. 

Ink = O-nk + Ctnk- (Bl) 
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(Recall that the electron dispersion is symmetric under this transformation.) Specifically, we can write ank = 
{\/2)[fnk + fn-k] and Unk = (l/2)[/nfc — fn-k]- Inserting this decomposition into the kinetic equation in the presence 
of a dc electric field in the y-direction (and without microwaves), we obtain the two equations 



dank _ gnfc - g.°fc 
dk Tin 



d(Jnk Oink 

Here, we have used the inequality Tin ^ r* and the fact that the disorder collision integral vanishes for the sym- 
metric part a„k- In addition, we have rewritten the collision integral for the antisymmetric part as {da/dt)dis = 
—Oink/'''*{enk)- Inserting the second of these equations into the first, we obtain 

f)2 

(ei^dc) (enfc)-^p- - ■ (d6) 

Note that this equation reproduces the estimate of the characteristic electric field E^^. This follows immediately from 
the fact that the characteristic scale of the k dependence is a/i\- 

For £'dc ^ -Edc therefore find cr„fc ~ it"^,. This is the starting point of the calculation leading to the expression 
(|37ll for the dark conductivity ayy. In the opposite limit E^c ^ ^dc write Unk = o'n + ^Unk, where ct„ is the 
average of ink over k. Note that ank which determines the current and hence the conductivity is directly related to 
Sank- Then, Eq. (|B3|1 shows that in magnitude 5(Jnk ^ (-^dc/-^dc)^o'5^fc- As a result, we expect that heating reduces 
the dark conductivity accor ding to Eq. gOl) in Sec. ITTTl 



APPENDIX C: EXPLICIT CALCULATION OF DISPLACEMENT PHOTOCURRENT jy 

In this appendix, we derive the displacement photocurrent in the y-direction more formally. In order to derive the 
quantum version of the meandering equipotential lines, we consider the Schrodinger equation, including the dc electric 
field in the y-direction, in LL representation 

{nk\Ha\n'k') - elJnn'Skk' - eEd,{-i) fnn'{k - k') (CI) 

with 

f f ^y^' (fc - fc')"'-"e-('=-'=')'/4i;^'-" f ^^4^) if n'>n 

fnn'{k-k')={ ^'""'^/2 " V 2 y - ^^^^ 

(fc-A;')"-"'e-('=-'=')'/4i::r"' (^^4^) if 7i>n' 

Neglecting LL mixing, we find the Schrodinger equation in the quasiclassical limit 

d 

(^nki^nk + eEdci-i)-gj^1pnk = £^nk (C3) 

with 

\n£) = Y,4,nk\nk). (C4) 

k 

This is readily solved and gives the quasiclassical meander states 

^nfe = V'Onexp|i^ dfc'^-^i^l (C5) 

To count the number of such states, we assume periodic boundary conditions in the x-direction, ^nk+L^/ej^ — i'nk, 
and thus (with / G Z) 

= (C6) 
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where £i is measured relative to the LL energy. As the energies £i fall into the range eE^cLy. the total mimber of 
states is L^Ly/2'Ki\, in agreement with the LL degeneracy. Requiring 1 = {n£i\n£i), we find the normalized meander 
states 



Wnk = {nk\n£i) = \l -j—S- exp \ i 

LixJ-ly 



ff" £ _ pO 1 



To verify that these states are indeed the meander states, we evaluate the expectation value of y = id/dk, and find 



k 



with £ni = u!c{n + 1/2) + £i, in agreement with the classical expectation. 

Following the same arguments as for the displacement photocurrent in the a;-direction, the current in the ?/-direction 
can now be expressed as 



_^.photoI = -^Y.Y. \<e.^n'S, - Vimfi^nl) - f{Sn'v)Wnl ' ^n'/' - Uj), (C9) 

^ nn' £i,£if 

where the transition matrix element is given by 

h^e^^n's.f = 2n\{n'£i,\T,\n£i)f. (CIO) 
By carrying out the summation over n, n' for Wc » T F we get 

^.photoi = e Ac^ ^ h%s,^^+i£, m£i - £v + Alo). (Cll) 

The relevant transition matrix element is 

\rNe,^N+is,f = 2^ (^^^' \{N£i\U\N£,) - {N + l£i\U\N + l£,)f. (C12) 

After disorder averaging, the matrix element becomes 

\{N£i\U\N£i,)-{N+l£i\U\N+l£v)\'' 

Inserting this into the expression for the current, we obtain 

Performing the energy sums yields 

E lE4S>f"^«''^^'^'''^'^l''^(^'-^''+^'^) = -^l / -e^^^'^''''*^°^(«^)+^+'«^('=+^)p. (C15) 
£. £. , " 6-E/(jc J 27r 

ti,ti/ k 

The fc-sum can be turned into an integral which in the limit E^^ — > can be evaluated in the stationary-phase 
approximation. This gives the result 



Note that E^^ drops out of this expression. This can be interpreted as follows. The length over which the electron 
can jump between meander lines is proportional to l/E^c- On the other hand, the electron density along the meander 
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line is proportional to -Edc- Thus, the overall probability to jump is independent of the dc electric field. In this way, 
we finally arrive at 

for the displacement photocurrent in the y-direction. Up to the sums over Landau levels, this is just the expression 
for the displacement photocurrent in the y-direction in Eq. (|55|) . 
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